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This  thMts  «xc«nds  th*  works  of  Noroo*  Laigheon,  and  Jonoo  idileh 
tako  a function  daflnad  by  a second  order  dlfforontlal  aquation  and 
datamlno  an  Interval  on  which  that  function  either  haa  a saro  or  attains 
a alnlata*  bounded  value.  Theoreas  for  locating  soroa  are  proved  for 
functions  of  a single  roal  variable,  functions  of  two  real  variables, 
and  real  vector*valuad  functlona.  Algorlthas  suitable  for  coaputer 
adaptation  are  presented  In  sasaplea  which  locate  seroa  of  such  fune- 
tiens  occurring  In  the  engineering  sciences  as  Legendre  polynoalals, 
laguarra  polynoalals,  toden's  aquation,  and  J)ufflng*s  aquation;  special 
alphas  la  Is  given  to  the  seroa  of  Basaal  ftmetlona  Jn(x).  The  aethods 
developed  In  thla  work  are  useful  In  optlalsatlon  theory;  they  also  can 
be  used  to  obtain  good  Initial  approalaatlons  of  seroa  for  starting 
Iterative  algorlthas,  such  as  the  Nawton-laphson  aethod,  which  give 
acre  aaact  saro  values. 


APPROXIMATIOM  OP  ZEROS  OP  FUNCTIONS 


ARISING  IN  THE  ENGINEERING  SCIENCES 


1h«  location  of  soros  of  polynonlals  and  other  functions  has 
always  been  a natter  of  great  Inportance  In  the  engineering  sciences 
because  these  seros  reveal  much  infomation  concerning  the  peculiarities 
of  the  problem  defined  by  the  finction.  Bessel  functions,  defined  by 
the  differential  equation 


y"  ♦ xy'  ♦ (x^  - p2)y  • 0 (y*-  ^ ) (l.l) 

QX 

have  a wide  range  of  applications,  from  their  first  use  by  Bernoulli  in 
determining  small  oscillations  in  a hanging  chain  to  the  modem  problem 
of  flux  distribution  In  a nuclear  reactor  (Ref  1x312).  Ihe  Legendre 
polynomials,  which  satisfy  the  differential  equation 


(1  • x^)y”  - 2xy'  ♦ n(nel)y  - 0 


(1.2) 


have  applications  In  problems  In  temperature  distribution  in  steadystate 

heat  flow  and  in  both  gravitational  and  electrostatic  potential,  as  well 

as  In  other  applications  of  a purely  mathematical  nature  (Ref  1x195). 

Another  application  of  finding  seros  of  functions  is  In  the  field 

of  optimisation  theory.  A basic  problem  in  optimisation  theory  is  that 

of  determining  the  variables  x^,  X2,  ...,  which  maximise  or  minimise 

a function  fOcj,  X2,  ...,  x^^)  (Ref  2x269).  To  do  this,  the  partial 

derivatives  are  set  equal  to  sero  and  the  resulting  n equations 

dxi 


1 


In  n taknoimn  am  solvad.  In  othar  words,  cha  alailtanaous  aaroa  of 

Xji)  “ (4  • l,2,»,.,n)  (1«3) 

am  to  bo  dotaralnod. 

Ihwro  sKlat  nany  ttaratlvo  sehsaMa  (tho  Nowton*&aphson  nathod  and 
tha  blaactlon  nathod,  for  axanpla)  which  can  ba  usad  to  find  tha  saros 
of  a function  f(x).  Howawar,  It  la  often  naeaaaary  to  start  tha  Itar- 
atlon  procasa  with  an  Initial  valua  idileh  la  alraady  eloaa  to  tha 
saro.  For  axanpla,  finding  tha  Imrarsa  of  tha  nuabar  2 la  aqulwalant 
to  solwing  tha  aquatlcn 

i . 2 • 0 (1.4) 

X 

ana  Itamtlwa  Mawten-Baphson  fotnula  for  solwlng  this  aquation  is 

Sal  • S^*  * 

An  Initial  valua  x^  • 2 uaad  In  Eq  (l.S)  producas  tha  dlwargant  sarlas 
I 2,  -4,  -40,  -3280,  ...  | ; Initial  valuaa  of  x^  - 1 and  Xo  - 0 pro- 
duea  tha  naanlnglaas  sarlas  1 1,0,0,...  | and  J0,0,0,...  | raspac- 
tlwaly.  Howawar,  If  Xq  • i la  uaad  to  start  tha  Itaratlcn,  tha 
rasultlng  sarlas  {1,  ...  J doas  convarga  to  tha  corract  solu- 

tlon,  X • y . 

Tha  problan  of  finding  good  starting  waluas  for  Itaratlwa  schsnas 
which  find  saros  of  a function  f(x)  la  consldarad  In  this  thasls  by 
datatnlnlng  an  Intorwal  [a,b]  whldi  contains  at  laaat  ona  saro  of  f(x). 
Than,  by  choosing  an  x^  fron  part»apa  at  tha  nldpolnt  or  an  end- 

point, tha  point  chossn  will  11a  "eloaa"  to  a saro  In  that  Intarwal. 
This  should  giro  a higher  probability  of  eonwargsnea  and  a faster  rata 
of  eonwargsnea  than  an  Initial  walua  Xq  chossn  at  randon. 


2 


trial  function  (and  tha  ona  uaad  pradoalnantly  in  this  chaptar)  la 


u(x)  • x(h>x)  • xh>x^ 


(2.2) 


Othar  adalsslbla  trial  functions  which  sight  ba  usad  on  £0»h3  ara 


slnl£ 

h 


(2.3) 


- xP  (h-x)** 


(2.6) 


vhara  P >!(<}>  1.  Corresponding  adnlsslbla  trial  functions  tsay  ba 
adapted  to  fit  any  Interval  [a,b]  or  [h,lch],  h > 0 and  k > 1,  to  extend 
the  dooMln  of  these  trial  functions. 

Associated  with  the  function  f(x),  the  trial  function  u(x)«  and  the 
differential  equation  L [f]  will  be  a functional  J £u3»  where 


J [**]  " / I A(x)  [u'(x)]^  - C(x)  [u(x;]^  I dx 


(2.5) 


Morse  (Ref  3 :32S>326)  proved  tha  following  theoreai  to  show  that,  when* 
ever  J ^u^  can  be  made  negative  by  varying  the  paraaeter  h,  then  a sero 
■ost  exist  on  tha  Interval  C^»h]]. 

Theores  2.1  * Morse  (Ref  3 ;325*326).  If  f(x)  Is  a function  which 
satisfies  Cq  (2.1),  and  If  u(x)  Is  an  adalsslbla  trial  function  on  Co,h] 
such  that  J [u]]  < 0 for  the  J [u]  defined  by  Cq  (2.5),  then  f(x)  has  at 
least  one  sero  on  [0,h]. 

Proof.  The  proof  by  contradiction  Is  begun  with  the  assumption 
that  f(x)  ^ 0 on  [|0,h].  Since  u(0)  > u(h)  ■ 0 by  definition.  It  can  be 


seen  that  the  following  Identity  holds: 


; r A(x)  f*(x)  u^( 
o 


] dx  - 


f - 0 (2.6) 

f(x)  O 


Also,  A(x)  > 0 and  ths  following  Insqusllty  holdst 


[u*(x)  ■ f > 0 


(2.7) 


Thsrsfors,  by  £qs  (2.6)  and  (2.7),  tha  following  Inoqualltlos  also  holdt 

0 < ; 1 [ y ^ . ujxl  {*(x)  ]2  I ^ 

“o'  f(x)  f 

or  (2.8) 

o</1[a(x)  f*(x)]’ aiiil  ♦ A(x)  f'(x)  [ iiiisi  ]’ 
o'  ^ f(x) 


♦ A(jc)  [u'(x)]*  - 2 A(x)  u*(x) 


♦ A(x)  [ afe)  f \ dx 

f (x)  ) 

which,  by  substitution  froo  Eq  (2.1),  Is  aqulvalent  to 
h 

I -C(x)  u^(x)  ♦ A(x)  f'(x)  C o (x)  ][ 

(x) 


0<  ; j -C(x)  u2(x)  ♦ A(x)  f’(x)  [ 2u,fr)  u*(x)  (,CxJ.  ■,  u2(x}  ] 

“ o ' f2(x) 


♦ A(x)  Cu»(x)]2  - 2a(x)  u’(x) 


♦ *(X,  j * 

r(x)  ’ 


Csncallatlon  of  taros  produces  tha  following  contradiction: 

h - 2 9 i 

0 < / I A(x)  [u*(x)  ] - C(x)  u^(x)  I dx  - J [u]  < 0 


This  contradiction  follows  froa  tha  assumption  that  f(x)  ^ 0 on  Co,h3. 
Thus,  f(x)  has  at  least  one  zero  on  tha  Interval  []0,h]. 


Exarole  2.1.  The  function  f(x)  ■ sin  x is  knovn  to  b«  • solution 
of  the  differential  equation 

♦ f(x)  - 0 (2.9) 


A conparlson  to  Eq  (2.1)  gives  A(x)  ■ C(x)  ■ 1.  Then,  with  the  trial 
function  u(x)  ■ xh*x^,  h Is  varied  until  J [u^  < Ot 

J [“D  " / I (h  • 2x)^  - (xh  • x^)^  I dx 


- 7 I (h^  • 4hx  ♦ 4x^)  - (h^x^  - 2hx^  ♦ x^)  | dx 


• (h^  • 2hx^  ♦ ^ x^  - ^ h^x^  ♦ j hx^  - y x*)  ] 


- h^  - 2h^  ♦ ^ h^  - j h^  ♦ h*  * 5 *'* 


- h’  < i - ^ h^X  0 


This  means  that  a sero  will  occur  on  [0,h^  when  ^ ^ ^ ^ 


h > /To  M 3.1623.  The  actual  sero  occurs  at  x > n <w  3.1616. 

Example  2.2.  The  function  f(x)  • Legendre  polynomial 

which  Is  a solution  to  the  differential  equation 


(1-x^)  f"(x)  -2x  f*(x)  ♦ n(n4l)  f(x)  - [(1-x^)  f'(x)]]  ♦ n(n4l)  f(x)  - 0 

(2.10) 

Now  A(x)  a 1 • x^  > 0 for  all  |x|  < 1,  and  C(x)  ■ n(n4l).  The  trial 
function  u(x)  > xh  - x^  and  Eq  (2.5)  are  used  to  determine  the  value  h 
for  which  J ^u]  < 0. 
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j C«]  • / [ • *^)0»  • n)(xh  - ] dx 


• / [ (h^-  Ahx  ♦ 4hx^-  4«^)-(n*#  n)(h^**^  2hs^*  xS  ] dx 

e 

- 0»^  - 2hx*#  I x3-  li^xal  Yi\^^  siaatl  hx*.  xh  f 

^3  2 S o 

• ^ - fel  t^Oin^o  n o 1)  - l5(n^o  n ♦ 2)  ♦ 60i^4  n ♦ 4)] 

3 30 

3 

■ flO  - h^  (n^o  a ♦ 4)3  < 0 
30 

Ihoroforo*  • »oro  will  occur  on  C®#**!  ***•“  ^ ( v h < 1, 

n^4  a o 4 

Idblo  I eoaparoo  tho  Malloat  pooltlvo  soro  of  P||(x)  to  tho  waluo  h for 
•onroral  intogor  valuoo  n > 3,  with  tho  true  soro  bolng  approxlaotod  f 
tho  gviorollsod  Rodrlguoo*  fonmlo  (Rof  4t8S2) 


?_(x)  - — i 1^2  (x*  - if 


2“  a!  -<“> 


Toblo  1 

Approclaatlon  of  Zoros  of  Logaidro  PolTnontols  P||(x) 


a.ii) 


Actual  Zoroa 


0.7746 

0.3400 

0.S38S 


ApprcKlaoto  Zoros 


0.7906 

0.6455 

0.5423 


0.3243 


0.3262 


It  !•  oft«n  of  Intoroot  to  find  soroo  of  f(x)  on  Intonrolo  not 
containing  the  origin.  Horo  generally,  the  Interval  h > 0 and 


k > 1,  nay  be  eonaldorod.  By  changing  the  Halts  of  Integration  and 
adjuatlng  the  trial  function  so  that  u<h)  • u(kh)  • 0,  the  following 
corollary  aay  be  established  In  the  ssae  aanner  as  Theorea  2.1. 

^oyo^lsQ.  If  f(x)  Is  a function  which  satisfies  Eq  (2.1),  and  If 
u(x)  is  sn  adkilsslble  trial  function  on  [h,kh]  su^  that  J [u^  < 0 for 

kh  9 , 

J M ■ ^ Cu'(x)]^  - C(jt)  u*<x)  J dx 

h 

then  f(x)  has  at  least  one  saro  on  Qh,kh3. 

In  this  case,  a convenient  trial  function  la 

u(x)  • (X  • h)(kh  - x)  - (k  o 1)  hx  - x^  - kh^  (2.12) 

Other  trial  functions  nay  also  be  constructed. 

Exsaole  2.3.  The  function  f(x)  ■ Bessel  function  of 

order  n of  the  first  kind,  which  Is  a solution  to  the  differential 
equation 

x2  f"(x)  ♦ X f*(x)  ♦ (x^  - n*)  fU)  - 0 (2.13) 

For  X > 0,  Sq  (2.13)  nay  bo  rewritten  to  natch  the  fom  of  Eq  (2.1)t 

9 2 

X f''(x)  ♦ f*(x)  ♦ (x  - £ ) f(x)  - [x  f'(x)]*  ♦ (x  • £-)  f(x)  - 0 

(2.10 


I 

.1 


Thus,  A(x)  • X and  C(x)  ■ x • ^ ; using  u(x)  froa  Cq  (2.12)*  J Is 
esleulstsdt 

kh  9 2 ««  <»  2 

C"]  ■ / I * C(M)  h - 2*3  * (x  • ^ )[(k#l)  hx  • - kh^3  | <** 


(2.15) 


- / { X [4x2  . 4 hx  ♦ (k4l)^h^] 

h 

2 

- (g  . 0.  )[x^-  2(k4l)hx2  ♦ (jk^*  4k*l)h2x2 

X 

- 2(k2  ♦ k)  h^x  ♦ k^h^3  | dx 

• ( x^  - ^ (k^l)  hx^  ♦ (k^  ♦ 2k  ♦ 1)  ^ x*  ♦ i.  (k^l)  hx^ 

- i 0(2  4 4k  ♦ 1)  h2x^  ♦ I 0(2  ♦ k)  hV  - i k2h^2 
4 3 2 

♦ n^  [J  x^  - I (k«l)hx2  ♦ (jkK  4k4l)h^x^ 

- 2 (k*4  k)h^x  ♦ k^^  In  X ])]** 


-(h^  (k^-1)  - I (k4l)  (k^-l)  ♦ i (k2+  2k*l)  h^(k2-l) 

- h*  (k*-l)  ♦ I (k4l)  h®  (k’-l)  - i (k244k^l)h*(k^-l) 
♦ § 0(*4k)  h*  (k^-l)  - J k2h*  (k2-l) 

♦ n2  [J  (k^-l)  - I (k^l)  (k^-l)  ♦ ^ (k2^  4k^l)h^(k2.l) 

- 20(*^  k)  h^(k-l)  ♦ k2h^  In  k ] ) < 0 


9 


Coablnlng  and  slapllfylng  tha  abova  axpraaalon  for  J [u]  glvaa 


h > 


12n*  In  k a (k^-l)  [(2an2)(k2al)  - (4  a flh2)kl 


)' 


(k«l)^(k2-l) 

f 

(2.16) 

This  Inaquallty  cannot  aaally  ba  aolvad  for  k If  h la  glvan  flrat. 
Chooaing  k flrat,  howavar,  aakaa  It  ralatlvaly  aaay  to  aolva  Eq  (2.16) 
for  h.  Ihua,  to  dotanlna  a apaclfle  for  an  Intanral  C^o'^o^o^* 
aaaaral  k aaluaa  mat  ba  ehosan  axparlaantally  until  tha  k^  la  found 
which  producaa  h^  In  Cq  (2.16). 

It  la  intoraatlng  to  axaalna  the  also  S(h,k)  of  tha  Intarval  ^h,kh3 
In  tha  Inatancaa  whara  k dacraaaaa  to  1 (and  h bacoma  vary  larga)  and 
whara  k bacoaMa  wary  larga  (and  h dacraaaaa  to  0).  Eq  (2.16)  aay  ba 
rawrlttan  to  axproaa  S(h,k),  tha  alaa  of  [h,kh^t 


S(h«k)  • kh-h  - (k'Dh  > ^3 


12n2k2ln  k ♦ (k2-l)C(2an2)(k2al).(4»8n2)k 


(k-l)2(k2.1) 

(2.17) 

Aa  k approachaa  1,  both  tha  nuaarator  and  danoalnator  of  tha  right  alda 
of  Iq  (2.17)  approach  aaro.  Applying  L'Hoapltal'a  Rula  thraa  tlma  to 
tha  Mipraaalon  within  tha  bracaa  In  Eq  (2.17)  producaa  tha  following 
Halt  for  S(h,k)t 


11a  S(h,k)  • 11a  (kh-h)  - /lo 

feol  kwl 

(h>b«»)  (h-w«») 


(2.18) 


Thua,  for  any  n,  aa  h gata  wary  larga  (and  k dacraaaaa  toward  1),  tha 
alaa  of  tha  Intarwal  containing  a aaro  of  f(x)  * J||(x)  will  bac< 
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elo««  to  JTo,  This  is  sxprsssod  aoro  sxsetly  by  tho  following  corollsryl 

Corollsrr.  For  sny  n > 0 snd  € > 0 thoro  oKlsts  a p > 0 sueh  thsc* 
for  sll  rosl  valuos  x > p,  thoro  Is  st  losst  ono  soro  of  JnCx)  on  tho 
latonrsl  [x,(x  ♦ •/  10)  ♦ €]. 

Ihoso  rosults  closoly  eorrospond  to  tho  fact  that  tho  dlstaneo  bo« 
twoon  soros  of  Jn(x)  approachos  n froa  bolow  for  n < % and  approachoa  n 
froa  abovo  for  n > % as  x grows  wary  largo  (Rof  S :49). 

Noxt,  by  allowing  k In  Eq  (2.17)  to  grow  Inflnltoly  largo,  tho  Halt 
of  tho  right-hand  sldo  again  bocoaos  flxod  for  a giwon  nt 

11a  S(h,k)  • 11a  (kh*h)  - [s  (2  o n^)]^  (2.19) 

k-s--  k-o«* 

(h-oO)  (h-»0) 

This  aoana  that,  as  k grows  vary  largo  snd  h approachos  0,  tho  saallost 
posltiwo  soro  of  Jn(x)  Has  in  tho  Intonral  f 0,(10  o 3n^)  o 6^  for  an 
arbitrarily  saall  €.  Tablos  II  and  III  list  sovoral  Intorvals  contain* 

Ing  soros  of  Jn(x)  obtainod  by  varying  tho  valuo  k in  Eq  (2.16)  and 
coapotlng  tho  corrospondlng  h.  For  coaparlson,  tho  first  10  soros  of 
J^(x)  for  sovoral  valuos  of  n aro  assoablod  In  tho  saao  tablos 
(Rofs  6-10). 

McCann  provos  that  Jn,p(x),  tho  p*^  soro  of  lowor 

bound 

^ ^ <*•*  (2.20) 

Rn  algorltha  for  calculating  an  uppor  bound  for  J„^p(x)  can  bo  dorlvod 
using  Eqs  (2.19)  and  (2.16): 

(1)  Caleulato  an  uppor  bound  k|h|  of  Jn^i(x)  using  olthor  Eq  (2.19) 
or  a k|  which  provldos  an  ovsn  saallor  bound; 
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TabU  II 

Coaparison  of  Intaxvala  Containing  Zaroa  to  Actual  Zaros 
of  Baaaol  Functions  Jn(x)»  npB0,I,2,3 

latarvali  * h«kh 


k 

0 

1 

e 

n • 1 

n - 2 

n - 3 

m 

0 - 3.163 

0 - 3.873 

0 - 5.478 

0 - 7.417 

so 

.064  - 3.227 

.077  - 3.882 

.107  - 3.390 

.144  - 7.234 

20 

.166  - 3.329 

.193  - 3.916 

.264  • 3.300 

.331  - 7.023 

3 

.790  - 3.933 

.862  - 4.313 

1.050  - 3.232 

1.304  - 6.321 

3 

1.381  . 4.744 

1.661  - 4.984 

1.881  - 3.644 

2.199  - 6.599 

2 

3.162  - 6.323 

3.232  - 6.466 

3.436  - 6.873 

3.750-  7.501 

1.3 

6.324  - 9.487 

6.373  • 9.363 

6.324  - 9.788 

6.767  - 10.131 

1.1 

31.622  - 34.783 

31.637  • 34.801 

31.680  - 34.849 

31.751  - 34.927 

1 

I 


Actual  Zaros  J.  -(x) 


r 

P-2 

P-3 

p-4 

P-5 

p-6 

P-7 

P-8 

p-9 

p-10 

i 

3.320 

8.634 

11.792 

14.931 

18.071 

21.212 

24.332 

27.49 

30.63 

■ 

7.016 

10.173 

13.324 

16.471 

19.616 

22.760 

23.90 

29.05 

32.19 

2 

8.417 

11.620 

14.796 

17.960 

21.117 

24.270 

27.42 

30.57 

33.71 

3 

6.380 

9.761 

13.013 

16.223 

19.409 

22.383 

25.75 

28.91 

32.06 

33.22 

Table  III 

Coaparlaon  of  Intarrala  Contalnlns  Zeros  to  Aetuol  Zeros 
of  Bessel  nmctlons  Jn^)»  n^»3t  10,30 

Interralt  ^h,kh3  • h-tdi 


n <■  4 

n • 3 

n ■ 10 

n m so 

0 • 9.487 

0 - 11.619 

0 - 22.384 

0 . 111.849 

SO 

.184  - 9.216 

.223  . 11.263 

.436  - 21.822 

2.139  - 107.933 

20 

.444  . 8.893 

.341  • 10.834 

1.044  - 20.887 

3.137  - 103.132 

n 

1.392  - 7.964 

1.900  . 9.304 

3.346  - 17.731 

17.301  . 86.310 

2.379  - 7.740 

2.999  - 8.998 

3.336  - 16.010 

23.332  - 76.399 

4.130  . 8.302 

4.614  - 9.230 

7.428  . 14.837 

33.733  - 67.312 

B 

7.092  - 10.639 

7.489  • 11.233 

10.216  . 13.323 

40.611  - 60.917 

i 

1.831  • 33.037 

31.979  - 33.177 

33.026  - 36.329 

37.169  - 62.887 

Actual  Zeros 


D 

P-1 

P-2 

P-3 

P-4 

1^5 

p-6 

P-7 

P-8 

P-9 

p-io| 

B 

111 

11.063 

14.373 

17.616 

20.827 

m 

27.20 

33.34 

36. 

70 

B 

|]Q] 

12.339 

13.700 

18.980 

22.218 

28.63 

38. 

16 

10 

14.476 

18.433 

22.047 

23.31 

28.89 

32.21 

33.30 

38.76 

42.00 

43. 

23 

so 

37.12 

62.81 

67.70 

72.19 

76.44 

80.31 

84.46 

88.32 

92.09 

93. 

80 

I 
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1 


I 


1 

I 


(2)  Choo««  another  SMillar  k2  and  caleulata  tha  eorraapondlng 
valua  h2  uaing  £q  (2.16).  MJuat  thla  lC2  and  raealeulata  h2 
oncll  li2  ia  graatar  than  but  eloaa  to  tha  uppar  botaMl  k|h|  of 
tha  prarloua  Intarval.  Ihua,  ttia  naw  Intarval  C*‘2*^2^2l 
contain  at  laaat  ona  saro  of  Jn(x)  and  will  ba  diajolnt  fron 
tha  pravlooa  Intarval. 

(3)  lapaat  atap  (2)  until  p diajolnt  Intarvala  hara  baan  eraatad. 

Slnea  aaeh  Intarval  eontalna  at  laaat  ona  saro  of  J||(x)»  tha 
uppar  Halt  of  tha  p^  Intarval  C**p**^**p3  *>•  upP*f 

bo«md  of  J„^p(x). 

Tha  abova  algorltha  will  ba  usad  to  caleulata  an 
uppar  bound  for  tha  fifth  saro  of  J3Cx>.  £q  (2.19)  nay  ba  uaad 

to  dataralna  that  at  laaat  ona  posltlva  saro  axlsts  In  tha  Intarval 
h ■ - [0,(10  ♦ 3n2)^  . [o,7.417].  Howavar,  using  k|  ■ 4.1 

In  Iq  (2.16)  produeaa  tha  Intarval  • [l. 583,6.494]  which  has  a such 
asuillar  uppar  bound  for  J3yi(x).  Naxt,  by  succasalvaly  axparinanting 
with  Iq  (2.16)  to  got  k|  sudi  that  h|  > tha  following  valuas 

of  k|  with  tha  asaoelatad  Intorvals  [h(,k|hj]  ara  obtalnad: 


**2 

“3 

k4 

•'5 


1.5234 

1.330 

1.244 

1.194 


t l2  • [6.499,  9.900] 

* I3  - [9.923,13.199] 

I I4  - [13.233,16.463] 

I I3  - [16.528,19.735] 


Each  of  tha  intorvals  13.5  contains  at  laaat  ona  saro.  Ihorafora,  tha 
uppar  Halt  of  I3,  19.735,  Is  an  uppar  bound  for  J3^3(x).  HoCann's  lowar 
bomd  in  Iq  (2.20)  for  J3^3(x)  Is  15.221.  Tharafora 


U 


I 


15.221  < J3^5(x)  < 19.735 


(2.21) 


The  exact  aero  from  Table  II,  *13,3 (x)  19.409,  is  much  closer  to  the 

upper  bound  than  the  lower.  However,  considerable  calculation  is  in« 

a 

volved  in  determining  this  upper  bound.  Even  greater  precision  may  be 

achieved  for  the  upper  bound  by  using  more  precise  k|'s  to  reduce  the 

distance  (h|  - ^i-l^ial^  between  the  intervals  I|  and  Since  the 

true  distance  between  zeros  for  J^Cx)  approaches  TT  3.14  for  large  x 

while  this  algorithm  produces  a distance  of,  at  best,  ^ 10  at  3.16,  there 

will  always  be  an  error  of  tit  least  no  matter  how  precise  ki  is 

3. 14  * 

chosen. 

Second  order  differential  equations  idiich  do  not  match  the  form  of 
Eq  (2.1)  may  often  be  converted  to  that  form  by  using  an  integrating 
factor.  For  the  differential  equation 


p(x)  f”(x)  ♦ q(x)  f’(x)  ♦ r(x)  f(x)  - 0,  p(x)  4 0 


(2.22) 


multiplying  both  sides  of  the  equation  by  the  integrating  factor 


I(x)  - 


I [q(x)/p(x)]dx 
p(x) 


(2.23) 


yields  the  equation 


/ [,(X)/P(X)>  ^ „ 

p(x) 

(2.24) 

Thus,  the  integrating  factor  can  be  used  to  convert  Eq  (2.22)  to  the 
self«adJoint  form,  Eq  (2.24),  from  which  expressions  for  A(x)  and  C(x) 
may  be  obtained  and  used  in  the  functional  J £u]  in  Eq  (2.5).  The 
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«xpr«ssions  A(x)  and  C(x)  may  be  too  complicated  to  Integrate  easily, 
however,  and  thus  be  too  unwieldy  to  use  effectively  in  J ^u]. 

2.2  Zeros  of  More  Complicated  Second  Order  Differential  Equations 


The  introduction  of  two  non-linear  functions  of  v(x),  a(w)  and  f(w). 

Into  Eq  (2.1)  produces  the  differential  equation 

»•  M - It  [a(x)  a [w(x)]  ^ ♦ C(x)  f [wCx)]  - 0 (2.25) 

whose  solutions  have  been  treated  in  some  detail  (Refs  12:377;  13:26-28). 

The  addition  of  another  function  B(x)  to  this  equation  produces  the 
following  more  general  abbreviated  form 

**  M - ^ [a(x)  a(w)  ~ ] ♦ 2 B(x)  ~ ♦ C(x)  f (w)  - 0 (2.26) 

Methods  to  establish  the  existence  of  real  seros  of  the  function  w(x) 
which  satisfies  Eq  (2.26)  are  similar  to  but  more  involved  than  those 
used  previously  in  this  chapter. 

Theorem  2.2.  If  the  following  conditions  are  tine: 

(1)  the  functions,  A(x),  a(v},  and  f(w),  are  differentiable  with 
respect  to  the  variables  x and  w,  where  A(x)  > 0 and  v(x) 
satisfies  Eq  (2.26),  and  B(x)  and  C(x)  are  continuous  on  the 
Interval 

(2)  for  an  admissible  trial  function  u(x)  on  [s,b],  there  exists 
an  associated  differentiable  function  G ([u^  such  that  G [u(a)l 
• G [u(b)3  « 0 and  G [u(x)]  > 0 for  a < x < b; 

(3)  f(w)  ^ 0 if  w(x)  4 0; 

Li 
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Ualag  • shorthand  fotn  dsrivsd  by  dropping  tho  varlsblss,  Cq  (2.30)  nay 
b«  raduoad  aa  followst 


■ / I 7 ^ Ca  w*D  ♦Aaw’^^y^a  A(u*)^  - A u*  * ® * 

♦ A W^**?  - 2B  u*  ♦ B ^ ft  **■  ♦ QC  I dx 

4f  * ' 

- / { f (-2B  »•-  Cf)  ♦ A a w*  ( ^ S **G  ) * A(ut)2.  AA  w*  /(M* 

a * ‘ 

, * <£  ^<»'-)".  . 2.  u.  g'  . . G^>  . q«  j *. 

- j Cu,o]  ♦ / ( f <-2B  *•)  ♦ U a »•  5JL  - A_a.*lg  u.*,) 

— a (.  l-»-.fa£*.2Lr.G  4 4.1^  g^te*}S  « B &AJ£l  I dx 

4f2  f > 

^ 2 

. J [..G]  . / ! <,2<,  - 4G  f»)  ♦ ^ (g  a - 2gS  1 dx  < 0 

a ' 4f^  ' ’ 


(2.31) 


This  contradiction  la  based  on  asataptlons  (4)  and  (S)  and  on  wCx)  4 0 
C**bDi  (A)  and  (5)  can  ba  satlsflad*  tha  aasunptlon  that  w(x)  4 0 
■oat  ba  falsa.  Thus,  w(x)  has  at  laast  ona  aaro  on  [a,b3. 

Iha  aiaibar  of  rastrlctlons  placed  on  tha  various  functions  In  this 
thaeraa  does  lialt  tha  range  of  applications.  Several  types  of  problsas 


dx 
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do  oKisc,  howovor,  which  fit  tho  fora  of  L [w]  In  Cq  (2.2S). 

Exoaolo  2.5.  Tho  Losuorro  pol3rnaaials  doflnod  by 

Mtlofy  tho  socond  ordor  dlfforontlol  aquation 

X I^(x)  ♦ (1-x)  I^U)  ♦ n l^U)  • 0 (2.33) 

This  any  bo  rovrltton 

[x  I^(x)]*  - X I^(x)  ♦ n l^(x)  - 0 (2.34) 

to  aoteh  Iq  (2.23),  with  A(x)  - x,  B(x)  - . ^ , C(x)  • n,  w(x)  > ^0c)» 

a(w)  ■ 1,  and  f(w)  ■ w(x)  • ^(x).  Q(x)  • ^ will  aaka  tho  aatrlx  M(x) 

la  Iq  (2.28)  posltlwo  aoal-dof Inlto.  Finally,  If  u(x)  • (xoa)(b-x)  and 

G rul  u^,  than  g(u)  > il£  • 2u  and  conditions  (1*4)  of  Ihooroa  2.2  ara 

du 

satisflod.  Tho  abowo  Inforaatlon  Is  than  substltutod  Into  Eq  (2.27)  to 
satisfy  condition  (S)t 

b - 

J Cu,c]  ■ / |x  (aob>2x)^  - 2 (-  y)(x-a)(b-x)(aol»2x) 

a ^ 

♦ (y  - a)  C (*-•)(*>-*) ]^  I dx  < 0 (2.33) 

To  find  an  Intorwal  containing  tho  first  posltlwo  soro  of  L„(x), 
tho  waluo  of  a In  Cq  (2.33)  la  aado  to  approach  soro.  It  can  bo  soon 
that,  as  a docroasos  toward  soro,  tho  waluo  of  tho  Intogrand  In  Cq  (2.33) 
owaluatod  at  a will  also  approach  soro.  Thoroforo,  with  a approaching 
soro  and  b ■ h,  tho  Halting  waluo  of  Cq  (2.33)  bocoaos 
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r 


J • / I x(h«2x)^  ♦ x(xhox^)(h-2x)  ♦ (J  -n)(xh-x^)*  | dx  (2.36) 

o 

- / { (h*x-4hx^4  4x^)  • (hV-  3hx3*  2x^) 

o 

«(!d  *3.  t xS  1 X*-  nh*x^*  2nhx^<-  nx*)  | dx 
4 2 4 * 

• (y  h^x*-  y hx^»  x^*  i y hx^*  y x*»  -jy  h^x*-  ^ hx* 

• 1 a h2«3^  n 4 n 5\  1*' 

♦ ^x-^hx^yhx-yx)  J 


l" 

o 


^ ^ w3/i  3 2-  ^6/1  1 lx 

-h(j-j»l)*  - J ♦ ^)  ♦ h <16  - 15  ♦ 


..  w3/  a ..  n Hv 

♦ h (-  ^ ♦ y • j) 


- h^(r  - 


^ - h ♦ ^ h^)  < 0 

60  240 


Sine*  h > 0,  J Qu,g3  < 0 lapll**  that 


[40  - 4(1  * 2n)  h * h^]  < 0 


(2.37) 


or 


(2  ♦ 6n)  - V (2  ♦ 4n)*  - 40  < h < (2  ♦ 4n)  ♦ V (2  * 4n)^-  40 


or 


h > (2  * 4n)  - V (2  * 4n)2-  40 


(2.38) 


(2.39) 


which  giv**  a r*al  valu*  h for  int*g*r  n > 2.  Tabl*  IV  coapar**  th* 
first  and  soeond  s*ros  of  s«raral  Laguorr*  poiynoaials  for  n > 2 with 
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th«  h froa  Cq  (2.39).  It  can  b«  accn  froa  the  table  that,  although  tha 
h occurs  wall  above  tha  first  saro.  It  does  fall  between  tha  first  and 
second  saros  In  tha  axaaplas  ahowi. 


Table  IV 

Zeros  of  Laguarra  Poljncalals 


n 

Polynoalal  ln(x) 

Interval 
with  Zero 

First 

Zero 

Second 

Zero 

2 

X x^  • 2x  V 1 

2 

[0,  2.254] 

0.586 

3.414 

3 

. 1 — x^-  3x  ♦ 1 

6 2 

[0,  l.5i0] 

0.416 

2.294 

4 

It  x^-  \ x^v  3x^-  4x  V 1 

24  3 

[O,  1.148] 

0.323 

1.746 

(Kef  4t8S3) 


2.3  Mlnlaua  Values  of  Certain  Second  Order  Differential  Equations 

In  addition  to  the  proof  of  sKlstence  of  a sero  on  an  Interval,  It 
aay  be  useful  to  deteralne  the  axlstanca  of  a alnlaua  value  of  a function 
v(x)  on  an  Interval  £a,b]],  where  v(x)  either  cosms  "eioae"  to  or  crosses 
the  X axis.  If  tha  Interval  [a,b^  contains  no  saros  of  v(x),  that  Is, 
either  v(x)  > 0 or  v(x)  < 0 for  e < x < b,  than  the  Interval  Ca,b^  la 
said  to  be  sere*free.  Kosdeov  derived  the  following  thoorsai  ccncomlng 
sero*free  Intervals  of  a function  v(x)  which  satlsflea  a certain  class 
ef  differential  equations. 

Theoren  2.3  - Konkov  (Ref  13t26-28).  If  a finctlon  v(x)  satisfies 
the  differential  equation 

L [v]  - [a(x)  v»(x)]'  ♦ C(x)  f Cv(x)]  • 0 (2.40) 
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•ttbjaet  to  tho  foilovlag  eonditlonot 

(1)  f C^U)]  - [v<x)]"  , n > 1; 

(2)  A(x)  > 0 Is  dlfforsntlablo  and  C(x)  la  continuous  on 

(3)  an  atelsaibla  trial  function  u(x)  on  [otb]  has  an  asaoclatod 

dlfforantlabla  function  C fu]  such  that  G ^u(a)3  ■ G ^u(b)3  ■ 0 
and  G ^ ^ a < X < b;  and 

(4)  J [u,g]  < 0 with 


b 2 

J Cu,g]  • / I A<x)  [u*(x)]  • CCx)  G [u3  J dx 


a.4i) 


than  ttiora  exlsti 


x6  [.,b]  * «L»J 


dG 

ilua  n ■ nax  _ i > 0,  g(u}  • such  that 


▼<X) 


. nlr 
< (§) 


<2.42) 


on  soaa  aublntorval  of  [a^bj. 

tJSSL*  ^2  v(x)  • 0 at  sons  point  In  Ca,b3»  than  tha  thaoran  la 
laasdlataly  trua.  Tharofora,  If  It  la  assunad  that  Ca,b3  la  saro«fraa, 
than  f (w)  • [v<x)3''*  ^ 0 on  [a^b]  and  tha  following  Inaquallty  may  ba 
darlrodt 


(2.43) 


f [(Ar»)*  G ♦ Aw*  gu*]  - Av*  G(|| 


-]) 


-(Cf)C  ♦ A KU*  A V*  G (nv""^)  v*  -,  ) 

♦ L j — j } ^ 

m / I [a(u*)*-  C g]  ♦ <*^  - AG  nr"*')  ♦ (u*g«gu*)  | dJC 

• J Qu^g]  ♦ / r ^^2^  (g^  - AG  nr"*')  1 
* ! Ar  J 

Sine*  J Cii,g3  < 0,  this  InpllM  that  (g^  • AG  nr"*')  > 0 for  Iq  (2. A3) 
to  b*  peoltlr*.  This  is  oquiralont  to  saying  that,  for  son*  x in 


C^ot)]  < 


jSsSi 


AG  [uOO 


n “ x€[[a,b3  AG 


or 


rU)  < 


C2.AA) 


■q 


r -i"*^ 

In  tho  spoclfie  caso  «h*r*  n is  an  odd  intagor,  L^0c)J  > 0 and 
(2.A2)  bacoMS 


|tu)|  < 5)^ 


(2.45) 


Ptanola  2.6.  EMan's  aquation,  which  occurs  in  astrophysics,  nay 
ba  written  as 


[x^r'U)]*  ♦ Cr(x)]"  - 0 


(2.46) 
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UhcrMs  the  full  solution  to  this  squstlon  Is  unknown,  Ihsorsn  2.3  asy 
bo  usod  to  dotonlno  how  eloso  ths  function  v(x)  will  cesio  to  soro  on 
£a,b3*  For  h > 0,  a ralua  p will  bo  usod  to  nlnlalso  h,  whoro 
G ^u]  • p u^,  u > x(h-x),  and  J []u,g]  < Ot 


jCu,g]  - / [a(u')*  - C g]  <toc  - / I x*(h^»  4hx  ♦ 4x^) 


- x^x^p  (h^-  2hx  ♦ x^)]]  } dx 


- / [ ♦ 4x^  - p(h^x^  - 2hx^  ox*)]  dx 


( . ^,4  . i,5  . ^,5  ^ . ^,7)  ]» 


* »’<  5 - 7 ♦ I ) • l>^<  I ■ f * ? ) 


This  lapllos  that  p > • and  tho  valuo  ■ bocoaos 


...  li  . . p 

x€  [0,h]  ^ x€  [0,h] 


(2.47) 


Thus,  In  tho  spoclfle  caso  whoro  n ■ S and  h • 2,  for  any  P > ^ ■ 3.5, 

2* 

w(x)  will  attain  a bound 


1 1 

\r  (x)l  < - (f)^ 


(2.48) 


an  tho  itttanral  [0,2j,  For  a p • 3.31,  tho  rosult  |v(x)|  < .916  Is 


slightly  bsttsr  than  Koakov's  raault  IvOc)!  < *943  for  ths  ssaa  «xaapls 
(Rof  12t380«381).  This  slightly  lovor  alnims  valus  of  v(x)  srlsos  froai 
tho  uso  of  s variable  coafflelant  p with  the  functional  G [u^. 

Another  ainlaua  value  theorea  aay  be  derived  froa  Iheoroa  2.2  by 
Introducing  an  intemedlata  tern  2B(x)  v*0c)  into  Eq  (2.40).  The  saae 
conditions  (1«S}  are  used  froa  Theorea  2.2,  except  that  a(w}  • 1 and  the 
expression  (g^(u)  - 4G  £u3  ) 1*  not  restricted  to  negative  values. 

Theorea  2.4.  If  the  function  v(x)  satisfies  the  differential 
aquation 

I.  [vj  • [a(x)  v*(x>]*  ♦ 2bCx)  v'(x)  ♦ C(x)  f [v(x)]  - 0 (2.49) 

and  If  conditions  (1-3)  of  Theorea  2.2  hold  on  an  interval  except 

that  g^(u)  - 40  £u]  ^ • K [u,f3  aay  be  poaitlve  and  that  a(v)  • 1, 
then  for  soae  x^€  [a,b]  either  (1)  f»  [vOtp)]  - — ^ < 1 or  (11)  v(x^) 

- 0.  * 

tsasL*  The  proof  parallels  that  of  Theorea  2.2;  assuaing  v(x)  4 0 
yields  the  following  Inequality! 

0 < J [u.c]  ♦ / r (g^  - 4f *C)  ] dX 

• L 4f^  J 

♦ ? [A-Sl-Si  (g  - 20^)  ] dx 

Since  (g  - 20^)  ■ 0 by  condition  (4), 

/ r (g2  - 4f  * G)  1 dx  < 0 

a L f2  J ” 


(2.30) 


iMda  CO  • contradiction;  this  iaplias  that  vOcg)  • 0 for  soaa 
Xg€  Howavar,  if 

b 2 

/ [ (g2  - 4f»  G)]  dx  > 0 

a ^ 

than  (g^  - 4f*  G)  > 0 for  sooa  Xg€  [a,b],  inplying  that  f • < ■ i. 

Thia  conplatas  tha  proof. 

Exanola  2.7.  A fora  of  Ouffing's  aquation, 

t"  ♦ p V*  a (q^y  ♦ rv^)  • 0 (2.51) 

which  arises  in  aachanics  in  tha  study  of  hard  springs  (r  > 0)  and 
soft  springs  (r  < 0),  nay  ba  invastigatad  using  Thaorsa  2.4  (Raf  14tl6-18). 
tfith  Q(x)  and  G • u^,  and  tha  trial  function 

u(x)  aXCh-x)  on  tha  intanral  [0,h],  J [u,G^  froa  Eq  (2.27)  bacoaaa 

J [u,g]  • / {a(u*)^  - 2B  G^  u*  ♦ (Q  - C)  G j dx  (2.52) 

o ' ' 

• / ^l(h  - 2x)^  - p(xh-x^)(h-2x)  ♦ ( ^ - l)(xh-x^)  j dx 

■ / [ h*  - 4hx  ♦ 4x^  - ph^x  ♦ 3phx^  - 2px^ 
o 

♦ ( ^ - 1)  (h*x^  - 2hx^  ♦ x^)  ] dx 

• - 2hx^  ^ 1 " 2 

♦ (£.  l)  (£.3  ]" 

4 3 2 5 -*0 
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- (h^  - 2h^  ♦ I h^)  ♦ (-  J ♦ ph^  - J h*) 


2 _2 


■ Cl  ♦ si  (4“  • 1)]  < 0 


This  larliM  that,  for  |p|  < 2,  h > 1 » tha  coofficlont  p 

datanilaoa  tha  sisa  of  tha  Intanral  undar  eons Idarat Ion.  If 

|p|  >2,  tha  diffarantlal  aquation  Eq  (2. SI)  nay  ba  muitlpliad  by  a 
constant  € to  naka  ]€p]  < 0.  This  will  changa  tha  coafficiants  A,  B, 
and  C in  L in  Eq  (2.49). 

Aftar  tha  sisa  of  tha  intarral  (|0,h3  is  dataminad,  tha  darivativa 


ft  (w)  - • 1 (q^T  a rr^)  - q^  ♦ 3 nr^ 

dr  dv 


(2.53) 


is  ananinad  to  datanina  if  v(x)  > 0 or  if  v(x)  conas  closa  to  saro. 
For  oxanpla,  if  p ■ 1,  than  h > (^)  sw  3.652  and  tha  intanral  undar 
SKSBinatian  is  £0,  3.652]].  If  q « r • j , than  7 + 7 v^Oc)  < 1 
inplias  that  w^(x)  < ~ , or  |v(x)  | < .7071  for  sobm  x€  [O,  3.652]. 
Howavar,  if  q « r ■ 1 in  Eq  (2.51),  than  it  is  always  tha  casa  that 
f*  • 1 ♦ T^(x)  > 1 sinea  v^(x)  > 0.  This  naans  that  - 4f*G  < 0 in 
Iq  (2.50),  which  inplias  that  wOc)  > 0 for  sons  x€  fO,  3.652].  Thus, 


Ihaoran  2.4  can  prorida  a boundad  valua,  and  in  sons  casas  a saro,  of  a 
fiaetion  w(x)  which  satisfies  tha  differential  aquation  L [w]  • 0 in 
Eq  (2.49). 


II 

L 
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1h«  abov*  th( 


la  also  applleabla  to  tha  aquation 


o"(t)  a (u  - ku^)  • 0 (2. 34) 

ahl^  is  usad  In  tha  titaory  of  aquatorlal  satallita  orbits  of  an  oblata 
spharold,  whom  u mpmsants  tha  variation  In  tha  radius  of  orbit  and 
t mpmsants  tha  angular  var labia  (Raf  14i23>24). 


mLMrnnDrri 


In  this  ohiM>ter  aufficlant  conditions  ars  astabliahad  for  tha 


cistanca  of  a raal  aaro  (or  a raal  ainiatai  valua)  of  a continuous  raal- 


valuad  function  f(X|^pc2)  with  two  raal  rariablas  within  a ragion  R.  Iha 


thaoraas  ara  praaentad  for  only  two  variablas;  howavar,  they  can  ba  ex» 


taided  to  apply  to  funoticns  of  savaral  variablas. 


Dia  raal  saroa  of  will  ba  considared  within  tha  ragion  R 

defined  by  tha  rectangle  R > { (X|«X2)i  ; for 


other  regions,  basic  transfonnations  nay  ba  applied  as  in  tha  previous 


chapter. 


Tha  function  f(x2,X2)  is  required  to  ba  a solution  of  tha  differen- 


tial aquation 


2 2 


W " ^ ^ 5x7J^*iJ^l»*2>  1^7  ^ * C(jti,X2)Cf(>Ci,X2)J 

i»l  J«1  * J 


(3.1) 


where  s is  an  integer,  s > 1,  C(x2,X2)  is  continuous  on  R,  and  tha 
alsnants  •|jO(i,X2)  ara  differentiable  on  R and  fora  a syaaetric,  posi- 


tive sani-daf inita  natrix 


A (X|,X2)  - 


•ll  <»1**2>  *12  <*1**2> 

•21  *22 


(3.2) 


Ateissibla  trial  functions  u(x2,X2)  will  bo  used  which  ara  diffarantiabla 


an  t,  sera  on  tha  boundary  of  R,  and  positive  in  tha  interior  of  R. 


1 
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TVo  such  trial  functions  for  ths  R dofinsd  sbovs  aro 
u0ci**2)  ■ *l  (h-Xj)  X2  (k-X2) 


and 


u(xi,X2) 


nxi  ^ nx2  , 

sin  ( ) sin  ( —=■  ) 

h k 


(3.3) 

(3.4) 


lhaaa  trial  functions  vanish  on  tho  boundary  of  R and  are  positive  with* 
in  R. 

Associated  with  f (X2»X2)  and  the  differential  equation  L [f^  > 0 
will  be  the  functional  J |^u]  which  is  defined,  in  its  most  basic  form, 
by 

J M - / / { » S 2 aij(xi,X2)  ( IJ-  ) ( ) 

» ial  J«1  * J 

- CCXjtXj)  u^'CXjfXj)  } (tXjdX2  (3.5) 

As  with  J [u^  in  the  single  variable  case,  %ihenever  J [u]  < 0,  there 
exists  at  least  one  sero  of  f(x2,X2)  within  the  closed  region  R for 
s ■ 1,  or  there  exists  under  certain  conditions  a minimum  value  of 
f(X|,X2)  within  R for  s > 1. 

The  following  lennw,  which  will  be  used  to  prove  theoresm  later  in 
this  diapter,  establishes  the  equality 

B [u.f]  - d(j[u])  (3.6) 

where  d(j([u3)  is  defined  to  be  the  integrand  of  jfu]]  in  Sq  (3.3)  and 
I^tt,f3  is  an  energy  functional  defined  below  on  a region  R in  which 
f(X|,X2)  4 0. 
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If,  within  a given  region  R,  f(xi,X2)  * function  which 
Mtlsflee  Eq  (3.1)  end  does  not  vanish,  u(X|,X2}  Is  an  achilsslble  trial 
function,  the  Integer  s > 1,  and  the  coefficients  a^jCxifXj)  fora  the 
positive  seni-deflnite  natrlx  shown  In  Eq  (3.2),  then  the  following 
identity  holds: 


d(jCu])  - ECu,f] 


(f) 


♦ s [l  - ( Y ) * 3 H (u) 


2 2 
♦ S E 

1-1  J-1 


r a 

^x^  «-  ^2s-l  ‘IJ 


¥-1 


dx 


(3.7) 


where  d(j(^u])  represents  the  Integrand  of  j[u]  In  Eq  (3.S)  and  H(*)  is 
defined  by 


H(*) 


2 2 
Z Z 

1-1  J-1 


•ij 


(3.8) 


Proof*  Using  the  fact  that  a|2  - a2|  In  the  s>anetrlc  matrix  A, 
the  expansion  of  £Cu,f3  and  cancellation  of  terms  produces  the  desired 
Identity: 


^ ] H(  f ) 


(s-l)(  J )^*  H(f)  ♦ sCl  - ( 7 )^*‘^]  H(u) 


♦ 


(3.9) 


I 


t 


- {2»-l)  >22  <2)  ’^)  } « { (g-l) 


.2*.df  2 2 2 


, <..„  .22<J)  ‘(I^)  ) . { . S. 


1-1  j-i  ‘J'axj'ax/ 


- a • 


2a.2 


2a  2 2 


-)  2 £ ^*1 4 
f2a-l  1^1  dxj  IJ  dxj 


2 2 


* * ^1 A 


2a  2 2 


* A <*u  llj>] 


’ * ii  f - "'• '] 
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j,  1^4 

♦ d(J  [u])  - d(j(|u]) 

Ihua  S [u,f]  • d(jCu]). 

If  • region  R can  bo  found  for  which  • trial  function  u(X|,X2.) 

■akoa  < 0,  than  f(X|^2)  will  olthar  vanish  sanowhoro  In  R or  attain 
SOM  alnlauB  value  on  R»  depending  on  the  value  of  a.  for  a ■ 1*  the 
Integration  of  Eq  (3.7)  gives 

- SJ  ■[»*.*]  dXjdXj  - /jj/  [f^  H(J)  } dx^dxj 

2 2 2 

♦ / / { E E Wii  ^1^)1  } <l*|dX2  (3,10) 

R 1-1  J-1  "1  ^ ^ ^*J 

The  result  of  the  first  Integration  Is  non-negative,  since  the  positive 
seal-definite  aatrlx  A nakes  H(*)  In  Eq  (3.8)  non-negative.  The  result 
of  the  second  Integration  Is  sero,  since  the  trial  function  u(X|,X2)  - 0 
on  the  boimdary  of  R.  This  leads  to  the  following  theorent 

Theoren  3.1.  If  f(x|,X2)  Is  a function  which  satisfies  the  differen- 
tial equation  1^0  - 0 in  Eq  (3.1)  with  s - 1,  and  If  for  soae  region  R 

2 2 ^ ^ 

<«]  -IJ  { ..J  - C . } d.,«j  < 0 (3.11) 

then  f(X|,X2)  - 0 In  soae  subregion  of  R. 


2s  2 

f2*-l  1-1 


2 2 
e { s E E 

1-1  J-1 
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groof.  If  It  Is  ssstiMd  that  f(X|^2)  ^ 0 on  R,  th«n  C [u^f]  Is 
walladafinsd  on  R.  Iha  eoablnation  of  Eqs  (3.10)  and  (3.11)  glvas  rlsa 
to  tha  contradiction 

0 1 d*idX2  - jCu]  < 0 

lharafora,  tha  assuaptlon  that  f(X|»X2)  ^ 0 on  R is  falsa,  and  tha  proof 
is  coaplata. 

Bcanola  3.1.  A Laplaca  aquation  of  tha  fora 

fjj,  ♦ f yy  ♦ f(x,y)  - 0 (3.12) 

■atchas  Eq  (3.1),  with  aj^^  ■ a22  " 1 «nd  m^2  * *21  * ^ making  tha 
■atrix  A positiva  sani*daf inita  and  with  0*1  and  a « 1.  Tha  trial 
ftinction  u(x,y)  • x(h-x)  y (k-y)  fits  a ractangular  ragion  R ■ 

{ (x,y)l  0<x<h,  0<y<k}.  To  datamina  valuas  for  h and  k such 
that  tha  ragion  R contains  a saro,  h and  k ara  axpandad  until  tha 
functional  jj^u^  < 0.  Uius 

jC«]  • / / { I ♦ I - I <«)^  ) 

Ic  ti  9 2 2 9 2 

■ / / { [(h-2x)(ky-y  )]  ♦ [(hx-x^)(k-2y)]  - [(hx-x^)(ky-y^)]  } dxdy 

o o 

• / { (ky-y^)^  / (h^  - 4hx  ♦ 4x^  • h^x^  ♦ 2hx^  - x^)  dx 
o o 

♦ (k«2y)*  / (h^x^  • 2hx^  ♦ x^)  dx  } dy 
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I 


i 


- / { (ky-y^)*  *3  * h 

« 3 2 


♦ (k.2y)" 


♦ 


I 

5 


- / (h’.2h5*|h^.lh’,ih’.ih») 

♦ (k-2y)^  (1  h’  - i h*  ♦ 1 h*)  } dy 

- 4 ^ h*)  / (k^y^  - 2ky^  ♦ y^)  dy 

o 

♦ (^  h3)  / (k2  - 4ky  ♦ 4y3)  dy 

■*  o 

* ^ ^ ♦ ‘sS  5 '■’ 


. nil 

900 


(lOh^ 


« lOk^)  < 0 


Ihanforc,  for  positlTo  h and  k valuoa,  (lOh^  - h^k^  ♦ lOk^)  < 0;  this 
lapliM  that  for  a flxod  valua  h 


(3.13) 


For  roal  solutions*  h aust  bo  grsator  than  J“lOi  also*  tha  aysaatry  of 
tha  functional  L ^f^  lapllas  that  k > ./To.  In  tha  Instanca  whara  h • k* 
a aquara  raglon  R rasults*  and  tha  solution  of  Sq  (3.13)  glyas 
h > 2i^T  M 4.472.  Thus*  avary  solution  to  tha  Laplaca  aquation  (3.11) 
aust  hara  a saro  on  R * { (x*y)t  0 < x < 4.472  * 0 < y < 4.472  } . 
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M«Kt»  Eq  (3.1)  is  Invsstlgstsd  Cor  Intsgsr  valuss  of  s > 1,  and  a 
■iBimai  valua  is  decaminad  for  a ragion  R. 

Thaorsa  3.2.  If  f(X|«X2)  is  a function  which  satisfies  L []f]  > 0 
in  Iq  (3.1)  with  the  intagar  s > 1,  and  lf»  for  soaa  ragion  R,  the 
functional  J ^u^  in  Eq  (3.5)  is  negative,  then 

Cor  SOM  point  €R. 

Proof.  If  t(Xi^X2)  • 0 sonawhara  in  R,  the  proof  is  trivial,  since 
the  trial  function  u(x^,X2)  > 0 in  the  interior  of  R.  If  f(X|,X2)  4 0 
in  R,  than  E Cu,f3  is  vall-dafinad  in  R,  and  integration  of  Eq  (3.7) 
produces  the  inequality 

0 > J [tt]  - I J E [u,f]  dX2dX2  - /g/  ■[  ^ dXjdX2 

♦ y (S-1)(H)**  H(f)  dXjdXj 

♦ ^ 1 - (j)**  *1  H(u)  <iX|<>K2 

2 2 2m 

♦ //[S  2 jI-  ( H_  ai,  )]  dX2dX2  (3.14) 

R 1-1  J-1  f^^  ^ 

Hie  first  and  second  integrals  in  Eq  (3.14)  are  non>nagatlva  since  the 
integrands  are  non-nagativa.  The  fourth  integral  is  saro  since  n - 0 on 
the  boundaries  of  R.  Thus,  for  J [u]  to  be  negative,  the  third  integral 
■oat  be  negatives  this  inplies  that  [l  • (^^**^]  < 0 in  som  subregion 
of  R. 
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Th«r«fore, 


.uv2s-2 


(X)‘-  ‘ > I 


£(Xj,X2)|  < max  u(X|,X2} 
(x,y)€R 


for  acme  point  (x2tX2}  €R,  and  the  proof  is  complete. 


A non-linear  Laplace  equation 


«XX  ♦ fyv  ♦ [«Oc.y)]  - 0 


(3.15) 


fits  the  differential  equation  L [f]  ■ 0 in  £q  (3.1)  with  ajj  - a22  - C « 1 
and  aj2  • «2i  " ® with  s • 2.  Eq  (3.5)  and  the  calculations  of 
Example  3.1  are  used  to  determine  an  R for  uhlch  J^u]  < Ot 

k It 

♦ h^k^)]  - / / (1  u2®)dxdy 


- ^ (h^^  ♦ h^k^)  - / / [(hx-x^)(ky-y^)]  dxdy 

^5  « rt 


k h 

■ tI'  (hV  ♦ h^k^)  - / (ky-y^)^  / (h^x^-4h^x^e6h^x®-4hx^4«®)dx]dy 


• ^ (h\^  ♦ h^k^)  - J(ky-y^)  { (^\^-|h^®<h|h^x^-  jhx®-*^®)]  } dy 


- Jr  (hV  ♦ h^k^)  - -il-  ; (k^y^  - 4kV  ♦ 6k^y®  - ^ky’  ♦ y®)  dy 

45  630  g 


■ fr  ♦ hV)  . ^ • -bid-  (8820h^  - h®k®  ♦ 8820k^)  < 0 

45  630  ^ 396900 


If  h > k,  a square  region  R results,  and  (8820h^  • h®k®e  8820k^)<0 
becomes  h^®  > 17640  or  h > (17640)^^^®  w 2.66.  In  this  square  region 
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IV.  Zto»  of  R—l  V«ctor-V«lu«<l  Functions 


In  optlaisatlon  th«ory,  • knovladg*  of  coros  of  the  derivatives  of 
veecor^valued  functions  is  needed  in  order  to  Baxiaise  or  ninlnise  those 
functions  subject  to  constraints.  In  this  chapter,  sufficient  conditions 
are  established  for  the  sKlstonce  of  a real  sero  of  a veetor*valued 
function 


(4.1) 


in  a f inite-diaenslonal  vector  space.  For  the  siopllf ication  of  notation 
the  letters  u,  v,  w,  and  y will  represent  vectors  u(t),  v(t),  w(t)  and 
y(t)  of  a single  real  variable  t,  A and  C will  represent  continuous  func- 
tions A(t)  and  C(t),  K will  represent  a constant  syssietric  positive- 
definite  n-by-n  suitrix,  and  other  letters  will  represent  constants. 

The  notion  of  an  inner  product  will  be  used  extensively  in  this 
chapter.  The  inner  product  of  two  real  n-dinensional  coluan  vectors 
V and  w shall  be  defined  by 


< v,w  > 


<(vj,V2,... ,Vn)  • (*itW2*** • »*n^  ^ 

v*w  - vjwi  ♦ V2»2  ♦ •••  ♦ 


(4.2) 


Osins  this  definition,  three  basic  properties  for  Inner  products  follow 
(lef  lSt6-7)t 
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(1)  < v.w  > - < ir,v  > i 

(2)  < CjV  ♦ C2W#  y > ■ cj  < v,y  > ♦ Cj  < w,y  > ; and 

(3)  < r,v  > > 0 If  V ^ 0,  and  < v,v  > • 0 only  If  v > 0. 


Tha  following  laanaa  aatabliah  basic  innar  product  idantitias  to 
ba  usad  latar. 

diffarantlabla  vactor-valuad  functions  v and  w, 
< v,w  >'  « ^ ( < v,w  >)  - < v,w'  > ♦ < V*  ,w  > (4,4) 

<1C 

££Oo£.  This  is  easily  established  by  using  tha  definition  of  an 
hmar  product  from  Eq  (4.2)  t 

< v*w  >*  - (^1*1  ♦ ^^2*2  ♦ •••  ♦ ^n*n^* 

- (ri*!*  ♦ vj’wj  ♦ V2W2’  ♦ ^2**2  ♦ — ♦ v„w„’  ♦ v„’w„) 

- (▼!*!*  ♦ V2W2'  ♦ ...  ♦ v„«ij") 

♦ (vi*wi  ♦ ▼2'''2  ♦ •••  ♦ Vn'*n> 

• < T.W*  > ♦ < V*  ,W  > 

If  R is  a constant  syomwtric  matrix  and  ▼ is  a vector 
fUBOtlon*  than  < v,Rv  >'  ■ 2 < y,Rv*  >, 

Proof.  Since  R is  constant  and  synutrie  and  tha  properties  in 
Iq  (4.3)  apply*  tha  following  identity  holds t 

< v*,Rv  > • (v*)*  (Rv)  - [(v')*r]v  - [r^(v')]^v  - [r(v')]*v 

• < tv*,v  > ■ < v,Rv'  > (4,5) 
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■ < ▼,av'>  ♦ < v.RT*  > 

• 2 < ▼,*»*>  (4.6) 


To  loe«eo  tho  soros  of  a voctor-valuod  function  v(t),  it  will  bo 
roqulrod  that  v aatlaflaa  tho  Inner  product  dlffarontlal  aquation  L 
of  tho  typo 

L [v]  - [a  < v.Rw  >*]'  ♦ C < v,Rv  > - 0 (4.7) 

wharo  A(t)  la  dlffarantlabla,  A(t)  > 0,  and  C(t)  is  continuous  on  an 
Intarval  Ca,b].  A vactor-waluad  adaisslblo  trial  function  u(t)  will  bo 
osad.  In  a oumnar  slailar  to  tha  trial  functlona  of  previous  chapters, 
with  tha  conditions  that  u(t)  Is  differentiable  on  Ca,b],  u(a)  > u(b)  > 0, 
and  u(t)  > 0 for  t€(a,b).  Associated  with  the  vector  u will  be  a func« 
tienal  J Qu^,  where 

b 

J C«]  - / [a(  < u,Ru  >’)^  - C(  < u.Ru  >)^]  dt  (4.8) 

a 

The  following  theorns  deaonstrates  that,  whenever  J [u]  In  Eq  (4.8) 
can  be  nade  negative  by  varying  the  Halts  of  Integration  a and  b,  v has 
at  least  one  sero  on  Cs,b]. 

Iheorea  4.1.  If  v(t)  Is  a vector  which  satisfies  Eq  (4.7),  and 
If  u(t)  is  a vector*valued  acbalsslble  trial  function  on 
J < 0 for  the  J ^u^  defined  In  Eq  (4.8),  than  v(t)  haa  at  leaat  one 
saro  on  Ca,b]. 


Proof.  Iho  proof  by  canCrodleClon  ia  dono  in  aueh  tho  amim  way  4 
la  tho  proof  of  Ihooraai  2.1.  Sineo  R la  a posltive>dof  Inita  ouitrix, 

< ▼•Rv  > > 0 by  dofinltion.  If  it  la  aaauoad  that  < v,Rv  > ^ 0 on 
Qa»b3»  than  tho  abovo  laanaa  nay  bo  uaod  to  produce  the  following 
contradiction: 


0 - / [( 


^ < u.Ru  A < v.Rv  >* 


< v.Rv  > 


)’]  dt 


(4.9) 


j j(  < u,Ru  > ^ < v,Rv  > ^ ^ < u^Ru  > ^ ^2  j 

^ ( < v.Rv  > ^ < v.Rv  > ) 


b 2 • 

• < v.Rv  > y 

a < v.Rv  > 


♦ A ( " ^ ^ > ♦ < u.Ru  >*  < v.Rv  > \2  ^ 

< v.Rv  > 

• / [“(  ^ ^ ^ ^ < v.Rv  >*  < u.Ru  > j2 

- < v.Rv  > < v.Rv  > 


• 2A  ( ^ u«Ru  > < u.Ru  > < v.Rv  > < v.Rv  > j 

< v,Rv  >2 


♦ A ( < u.Ru  >*)  ] dt 


. j r/  < u.Ru  A < v.Rv  >*  1 . f < 

, < v.Rv  > / « ' < 

. 2A  < <“!«“><  >'  < v.«’'  >' 

< v.Rv  > 

w 2 ^ ^ 

• / r(  ^ ^ ^ ^ ^ )'  ♦ A ( > )2 

^ < v.Rv  > < v.Rv  > 


. A ( » < “i«V  > < "■>>“  >’  5 Ti^  >'  ) , A ( < .,««  >')^]  dt 
< v.Rv  > 


dt 


AS 


. J [tt]  < 0 

Hum*  « contradiction  arises  froa  the  asstaq>tlQn  that  < v,Rv  > ^ 0. 

This  asaas  that  < > • 0 scaoahara  on  [a,b3»  vhlch  lapllas  that 

v(t)  ■ 0 for  seas  t€  Qa^b^. 

Specific  exaaples  of  v(t)  which  satisfy  L In  Eq  (4.7)  are  dlffl> 
eult  to  find.  However,  differential  equations  already  In  the  fora  of 
■q  (4.7)  aay  bo  exaalned. 
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Ex— Pl«  4,1.  In  th«  CM«  whar*  v(c)  la  « scalar  function  of  t, 

E ■ 1,  tt(t)  • t(h»t)  on  the  Intanral  C^**']*  • 1 * and  C(t)  ■ 

n(B«l)  arc  eoafflclants  froa  tha  Lagandra  polynoalal  axaalnad  In 
Quvtar  II,  Eq  (4.7)  baconas 

L [r]  • [(1-t^)  < v(t),  vCt)  >']  ♦ n(nal)  < v(t),  v(t)  > - 0 (4.10) 

In  ordar  for  v(t)  to  hava  a saro  on  Co,h],  J [u]  daflnad  In  Eq  (4.8) 

■use  ba  nagatlva.  Thus 

h t 2 2 

J C**D  • / C(l-t^)(  < u,u  > ) - n(nal)  < u,u  > ]dt 
o 

• / < u,u'  >)^  » n(n4l)  < u,u  3dt 

a 

■ / C<l«t^)  4(uu')^  - n(nal)  (u^)^  ]dt 
o 

- / C(l-t^)  4(ht  - t^)^  - 2t)^  - (n^a  n)(ht  - t^)^  ]dt 
o 

. / [(4-4t^)(h^t^-  6h^t^a  I3h^t^-  I2ht^a  4t*) 
o 

- (h^a  n)(hS^-  4h^t*a  6h^t*-  4ht^a  t®)]dt 

- / [(4hS^-  24h^t^a  52h^t^-  48ht^a  I6t®) 
o 

- (4hS^-  24h^t^a  32h^t®-  48ht^a  I6t®) 

- (n^a  n)(h^t^-  4h^t®a  6h*t®-  4ht\  t®)]dt 
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- h*(n**  n)(y  - 

2 

" ^ 313“  ♦ *^55zr  " 6^  n ♦ 4)]  < 0 

Ih^rwfor*,  a saro  of  v(t)  will  occur  on  fO-hl  whan  h > (-»' 

n^w  n ^ 4 

and  h < !•  Whan  n ■ 3»  a saro  lias  In  ^0,  .8663;  whan  n • 4,  a saro 
lias  in  [O,  .707]. 

It  nay  ba  difficult  to  find  vector  functions  which  satisfy 
L ^v]  • 0 in  Eq  (4.7).  Howavar*  it  is  possible  to  substitute  the 
differential  inequality  L [v]  < 0 into  Eq  (4.7)  and  still  derive 
lhaoran  4.1  without  altering  the  proof  substantially.  This  allows  a 
wider  range  of  vector-valued  functions  v(x)  to  satisfy  the  resulting 
inequality. 

Other  results  by  Jones  have  been  established  using  different 
conditions  on  the  vector-valued  function  v(t)  and  the  vector-valued 
trial  functions  u(t)  (Ref  16:86). 
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V.  Conclusion  { 

Ihoro  aro  savoral  conputor  applications  of  itaratlva  tachnlquas 
which  will  Iterata  toward  saros  of  a function,  but  tha  initial  starting 
▼alua  for  thasa  schssias  oftan  nust  ba  chosan  carafully.  If  tha  Initial 
point  la  not  closa  to  an  actual  saro,  tha  Itaratlva  tachnlquas  often 
dlvarga. 

Tha  thaorsas  and  axaoplas  prasantad  In  this  work  can  ba  usad  to 
locate  Intervals  on  which  a function  has  a saro.  If  tha  given  function  | 

f can  ba  paired  with  one  of  tha  differential  aquations  L [f]  - 0,  and 
If  a trial  function  u Is  found  which  aakas  tha  corresponding  functional 
J [u]  nagatlva  on  a certain  interval,  than  a saro  of  f does  exist  on 
that  Interval. 

txanplas  of  cosMon  differential  aquations  tram  different  anglnaar* 

Ing  sclancas  ware  chosan  to  denonstrata  tha  breadth  of  applications 
possible.  Perhaps  this  nathod  can  ba  applied  to  Indicate  tha  existence 
of  saros  of  other  differential  equations  before  valuable  coaq>utar  tine 
la  wasted  In  hlt«or-nlss  saarchas  to  find  thasa  saros.  In  fact,  tha 
algorithns  obtained  can  thansalvas  ba  Inplanantad  using  nunarlcal  tach- 
nlquas an  a cooputar  of  reasonable  sisa.  ^ 


I 
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